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ABSTRACT 

We present the resuhs of a molecular line survey of the extreme carbon star 
CRL 3068. The observations were carried out with the Arizona Radio Obser- 
vatory (ARO) 12 m telescope and the Heinrich Hertz Submillimeter Telescope 
(SMT) at the A 2 mm and A 1.3 mm atmospheric windows. The observations 
cover the frequency bands from 130-162 GHz and 219.5-267.5 GHz. The typi- 
cal sensitivities achieved are Tr < 15 mK and Tr < 7mK for the ARO 12 m 
and SMT, respectively. Seventy two individual emission features belonging to 23 
molecular species and isotopologues were detected. Only three faint lines remain 
unidentified. The species C-C3H, CH3CN, SiC2, and the isotopologues, C^''0 and 
C^^O, HC^^N, HN^^C, C^^S, C^^S, ^^CS, ^sSiS, and ^^SiS are detected in this 
object for the first time. Rotational diagram analysis is carried out to deter- 
mine the column densities and excitation temperatures. The isotopic ratios of 
the elements C, N, O, S, and Si have also been estimated. The results are con- 
sistent with stellar CNO processing and suggest that CRL 3068 is more carbon 
rich than IRC+10216 and CIT6. It is also shown that the chemical composition 
in CRL 3068 is somewhat different from that in IRC+10216 with a more exten- 
sive synthesis of cyclic and long-chain molecules in CRL 3068. The results will 
provide valuable clues for better understanding circumstellar chemistry. 

Subject headings: ISM: molecules — radio lines: stars — line: identification — 
stars: AGB and post-AGB — stars: circumstellar matter — stars: individual 
(CRL 3068) — surveys 



1. Introduction 

As a star evolves up the asymptotic giant branch (AGB), copious amount of material 
at its outer layers are ejected into space through stellar winds and form an expanding cir- 
cumstellar envelope. Through millimeter-wave and infrared spectroscopy, we have learned 
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that circumstellar envelopes around evolved stars are very efficient chemistry factories that 
produce organic and inorganic molecules over very short time scales. However, our under- 
standing of the interactions between molecular processes, chemical environment, and local 
physical conditions is far from complete. In oxygen- and carbon-rich environments, the re- 
action routes are clearly different. In oxygen-rich objects, the circumstellar chemistry is 
believed to be dominated by OH and H2O molecules, and the che r nical complexity of the 
oxygen-rich environments has been demonstrated by IZiurys et al.l (120071 ). In carbon-rich 
objects, a large fraction of oxygen atoms is locked in CO and C-bearing compounds such as 
hydrocarbons and cyanopolyynes are efficiently reprocessed. Among the ca rbon-rich objects. 



the nearby AGB sta r IRC-l-10216 has been extensively studied (see, e.g. ICernicharo et al. 



2OOOI : iHe et al.ll2008l . and references therein), with more than 60 molecular species detected. 



This paper is one of a series of studies devoted to the understanding of circumstellar 
chemistry of evolv ed stars. In p revious work , we have presente d the spectra of the AGB 
stars IRC+10216 JHe et allboosh a nd GIT 6 Jzhang et allboooh and the planetary nebula 
(PN) NGG 7027 Jzhang et al.lboosh obtained using the Arizona Radio Observatory (ARO) 
12 m telescope and the Heinrich Hertz Submillimeter Telescope (SMT). In the present study, 
we report a molecular line survey of the carbon-rich circumstellar envelope GRL 3068 using 
the same telescope settings with the aiming at exploring the molecular constituents and 
chemistry processes in an extremely carbon-rich environment. 

The carbon-rich AGB star GRL 3068 (AFGL 3068, IRAS 2 3166+1655, LLPeg) i_s a long- 
period variable with a very long pul sation period of 696 days (ILe Bertre et al.lll995l ). It has 



a high G/0 abundance r atio of 1.38 (IWinters et al 



1997) and has been taken as a prototype 



of extreme carbon stars (jVolk et al.lll992l ). iJones et al.l (jl978l ) investigated the infrared (IR) 
spectrum of GRL 3068 and detected an SiG absorption feat u re at 11.3 fim. The presence 
of SiG absorption was further conffimed by iGlement et al.l (120031 ) based on the Infrared 
Space Observatory (ISO) observations. The fact that the SiG feature is in self absorption 
suggests that the circumstellar envelope of GRL 3068 is extremely optically thick. This is 
sup ported by the very low c o lor temperature o f the dust ena i ssion co mponent 300 K, see 
e.g. iLebofskv fc Riekel Il977l : IVolk et all Il992l : lOmont et all Il993h . iLe Bertre et^ Jl995h 



found that the dust continuum spectrum of GRL 3068 the wavelength range 1-50 /im is 
consistent with a dust opacity law of oc A~^ '^. 



Based on a consistent time-dependent hydrodynamical model, IWinters et al.l (119971 ) 
determined the distance to GRL 3068 to be 1.2 kpc, the time averaged outflow velocity as 
14.7 km s~^, and the average mass loss rate as 1.2 x lO~^M0yr^^. The distanc e derived by 



Winters et al.l (119971) is c l ose to a more recent value of 1.05 kpc reported by lYuasa et al. 



(I1999I ) and iMenzies et al.l (120061 ). The derived mass loss rate of IWinters et al.l (119971 ) is 2-10 
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times higher than t hose obtained using CO lines (IKnapp &: Morria Il985l : IVolk et aLlll993 



Schoier et a 



12002 



Woods et al. 2003 



Tevssier et ali boooh and tha t by fitting the dust 



contin uum (jVolk et aLlll992l ). The H I fine at 21cm was detected by iGerard &: Le Bertre 
( 120061 ). who est imated that the H I emis sion in CRL3068 is extended over a region 



8' in diameter. iMauron fc Huggind (120061 ) imaged the extended circumstellar envelope of 
CRL 3068 in the optical and were able to trace the galactic scattered light out to a distance 
of ~ 40" from the star. They estimated the expansion timescale of the envelope to be 
15,200 yr. The optical image was characterized by a remarkable spiral pattern traced out 
to 12" from the center, suggesting that binary companions might play a role in shaping the 
circumstellar envelope. The scenario of binary-induced outflow was further comflrmed by 
the near-IR adaptive optics images obtained with the Keck ll telescope, which showed that 
the ce ntral region of CR L 3068 is composed of two components separated by a distance of 
0.11" ( iMorris et al.ll2006l ). iNeri et al.l (119981 ) mapped the CO(l-O) emission, which revealed 
a compact inner envelope and a surrounding detached bipolar shell. 

Molecular lines in CRL 3068 at millimeter wavelengths have been ex plored by several 
autho rs. Relatively strong CO (1-0) line emission has been detec ted bv iKnapp fc Morris 
(Il985|) andhvman et all (1992h. HC3N lines were flrst detected bv ljewell fc Snvdeil h9m . 
Nguven-Ouang-Rieu et al.l (119881 ) discovered the HCN(l-O) and H^^CN(l-O) line emission. 
Sopka et all Jl989h reported t he observations of CO, ^^CO, HCN, H^^qn, HC3N, and HNC 
in this object. IWannier et al.l (jl99lh discovered relatively faint emission from C3H2. Radio 
observations toward CRL 3068 in the frequency ranges bet ween 39-47 GHz and 8 5-91GHz 
were carried out with the Nobeyama 45 m radio telescope by lFukasaku et al.l (Il994j ). and the 
molecular sp e cies C 2H, HC5N, C4H, and SiS were detected in this object for the flrst time. 
Woods et al.l (120031 ) reported a molecular line survey made with the SEST 15 m and Onsala 
20 m telescopes. The molecular species in CRL 3068 that were detected positively in their 
survey observation include CO, ^^CO, CN, HCN, H^^qn, HNC, C2H, HC3N, C3N, CS, and 
SiS. 

In this paper, we present spectra of CRL 3068, representing the most complete molec- 
ular line survey of this carbon-rich envelope so far. The remaining parts of this paper are 
organized as follows: the observations and data reduction are described in Sect. 2; in Sect. 3, 
we presents the identiflcation and measurements of detected molecular lines; Sect. 4 gives 
the column densities and abundances of detected species; in Sect. 5, the implication of our 
findings on circumstellar chemistry is discussed; our conclusions are summarized in Sect. 6. 



-4- 



2. Observations and data reduction 

2.1. The ARO 12 m observations 

The spectral hne survey at the A 2 mm window was carried out during the period from 
2007 December to 2008 January with the ARO 12 m telescope at Kitt Peak. The 2-mm dual- 
channel SIS receivers were employed, operated in single sideband dual polarization mode. 
The image rejection ratio is larger than 18 dB on most of the occasions. The system noise 
temperatures were typically 250-550 K. The spectrometer backends employed are two 256- 
channel filter banks (FBs) with a channel width of 1 MHz and a millimeter autocorrelator 
(MAC) with 3072 channels and 195 kHz per channel. The observations were made in beam 
switching mode with an azimuth beam throw of 2'. We checked the pointing every 2-3 hours 
by observations of planets. The pointing accuracy was better than 15" at normal weather 
conditions. The scanned frequency range was from 130.0 to 163.9 GHz with a few gaps. 
The half power beam width (HPBW) is about 40" at this frequency range. Typical on- 
source integration times ranged from 70-90 minutes for each frequency setting. For A 2 mm 
observations, w e gave priority to the fr equency ranges covering the strong lines as observed 



in IRC-l-10216 (jCernicharo et al.ll2000l ). Only a few extremely faint lines in the spectra of 
IRC+10216 fall within the gaps of our spectra. The temperature scale was given in term 
of T^, which was calibrated by the chopper-wheel method, and corrected for atmospheric 
attenuation, radiative loss, and rearward and forward scattering and spillover. The main 
beam brightness temperature was determined using Tr = T^/rj^, where t]^ is the corrected 
beam efficiency (~ 0.75). 

We used the CLASS software package in GILDAS for the data reduction. After 
discarding the bad scans, we obtained the co-adding spectra from individual scans. On some 
occasions, the MAC spectra might suffer severe effects from the bandpass irregularities, and 
thus for some frequency coverages the signal-to-noise ratios might be relatively low due to 
the exclusion of these scans. The baseline was determined by low-order polynomial fits to 
the line-free regions of the spectra. In order to improve the signal-to-noise ratio, the spectra 
were smoothed and rebinned to a frequency resolution of 1 MHz, corresponding a velocity 
resolution of ~ 2kms~^. The MAC spectra are presented in Figure [H The typical rms noise 
temperature is less than 15 mK in main beam temperature unit. 



^GILDAS is developed and distributed by the Observatorie de Grenoble and IRAM. 
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2.2. The SMT observations 

The spectral line survey at the A 1.3 mm window was made using the SMT 10 m telescope 
located on Mt. Graham in 2008 January. For the observations, we utilized the ALMA 
Band 6 (211-275 GHz) sideband-separating SIS mixing-preamp, yielding a typical image 
rejection ratio of > 20 dB and a system noise temperature of ~ 150 K. The 2048-channel 
acousto-optical spectrometer (AOS) and the 1024-channel Forbes Filterbanks (FFBs) were 
used simultaneously and were configured to provide a channel width of 500 kHz and 1 MHz, 
respectively. A beam switching mode was applied with a separation between the ON and 
OFF beams of 2' in azimuth. The pointing was checked every 2 hours on nearby planets 
and was found to be accurate to within 10". The complete frequency range from 219.5 
to 267.5 GHz was covered, with typical on-source integration time of 30 minutes at each 
frequency band. The beam size at this frequency range is about 30". The data were calibrated 
to the temperature scale, which was corrected for atmospheric attenuation. The main 
beam temperature was derived using Tr = T^/rimb, where rjmb is the main beam efficiency 
(~ 0.7). 

The procedure of data reduction was same as that for the 12 m data. Since there is no 
gap in the A 1.3 mm survey region, all the spectral bands were merged together. Figure [2] 
gives the FFB spectra, which have been smoothed and rebinned by a factor of 3, yielding 
a frequency resolution of 3 MHz, corresponding a velocity resolution of 3.7kms~^, and a 
typical rms temperature of 5mK in main beam temperature unit. 

3. Observational results 

3.1. The overall survey 

The spectra together with line identifications are shown in Figures [1] and [2l Note that 
some features in these figures are not real but caused by bandpass irregularities or bad chan- 
nels. A total of 75 distinct emission features (including a few uncertain ones) were detected 
in our survey observations. The most intense line detected in our survey is the CO (2-1) 
transition, followed by the HCN (3-2) transition. Apart from these two transitions, the 
others have integrated intensities of less than lOKkms"^. The molecular lines detected in 
the high frequency window are a factor of about two more than those detected in the low 
frequency window. This is partially due to the fact that the ARO 12 m telescope has a 
larger beam size than the SMT and thus the 12 m data are more likely to be affected by the 
beam dilution. The molecular line databases that are used for line identifications include 
the archives of molecular line frequencies derived from the theoretical calculations of the 
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JPL cat alog (Pickett et al 



rCD MS. iMuUer et al 



Cernicharo et al 



2001 



,1998l|^ and and the Cologne database for molecular spectroscopy 
l2005llN. and those from previous observations of othe r 



2000 



sources 



Pardo et al.l 120071 : iHe et al.l l2008l : IZhang et al.l l2009l . and the 



(e.g. 

NIST Recommended Rest Frequencies for Observed Interstellar Molecular Microwave Tran- 
sition^. We identified 11 main molecular species, including CO, CS, C2H, C4H, C-C3H, 
CH3CN, CN, HCN, SiCs, SiO, and SiS, and 12 rare isotopologues, including ^^CO, C^^O, 
C^^O, C^^^S, C^^S, i^CS, HC^^N, H^^CN, HN^^C, HC3N, ^QSiS, and ^^SiS. Only three fairly 
weak features remain unidentified. Several transitions are new detection in CRL 3068. A full 
list of line identifications and measurements is presented in Table [TJ Columns 1-3 give the 
identified species, transitions, and frequencies, respectively. The rms noise levels are given 
in column 4. Column 5 and 6 lists the peak and integrated intensities. Column 7 gives the 
line widths (FWHM) obtained by fitting Gaussian line profiles. For blended lines. Table [T] 
gives the intensities and widths of the combined features. 

Figure [3] shows the cumulative number of detected lines exceeding a given integrated 
intensity. For comparison, the figure also plots the results for the spectra o f IRC-l-10216 



and C IT6 which were obtained in the same observation project (IHe et al.ll2008l : IZhang et al. 



2OO9I ). The shapes of these curves in Figure [3] are similar with each other although that of 
CRL 3068 seems to be slightly steeper than the other two. This suggests that the differences 
between the three objects can be mostly a distance effect, and the detection of a large number 
of lines in IRC+10216 does not imply that IRC+10216 is unique or richer in molecular content 
than the other two. 



3.2. Individual molecules 



The details of the molecules detection in CRL 3068 are described below. Here we com- 
pare thes e dete ction s with the spectra o f IRC+10216 and CIT6, which are respectively from 
He et al.l (120081 ) and IZhang et al.l (120091 ) if no statement is issued. 



^http:/ /spec.jpl. nasa.gov. 
http://www.phl. uni-koeln.de/vorhersagen/. 
''http: // physics .nist .gov / cgi-bin/micro /table5 /start .pi 
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3.2.1. CO 



The J = 2-1 transitions of CO and its isotopologues ^'^CO, C^''0, and C^^O were de- 
tected by the SMT. The CO (2-1) transition is the strongest hne detected i n this survey, with 
an integrated intensity of 35.9Kkms~^. This hne has been observed by IVolk et al.l (119931 ) 



using the 15 m JCMT. They obtained an integrated intensity of j T\(lv = 40.2Kkms~^ (i.e. 
/ T^dv = 57.4Kkms~^), shght ly larger than our re s ult. A higher integrated intensity of the 
CO (2-1) line was detected by iGroenewegen et al.l (119961 ) using the IRAM telescope. How- 
ever, considering the larger beam size of the SMT compared to those of JCMT and IRAM, 
the se intensity discrepancies are i nsignificant. The ; ^ ^C0 ( 2-1) line has also been detected 
by IGroenewegen et al.l (119961 ) and iBujarrabal et al.l (Il994j ) . Our result yields an intensity 
ratio of /(^^CO 2-l)/I{^^C0 2-1)= 0.14, larger than the values of 0.10 and 0.07 found in 
IRC+10216 and CIT6. This reflects differnt ^^C/^^C isotopic ratios and/or different optical 
depth effects on the three objects. The faint C^''0 and C^^O lines are only marginally above 
the detection limit. The two rare isotopologues should be new detections for this object. 
The CO line shows a narrower profile than the isotopic lines which can be explained by the 
fact that the former is optically thick. 



3.2.2. CS 

The J = 3-2 and J = 5-4 transitions of CS were detected by the ARO 12 m and 
the SMT, respec t ively. The two CS lines toward this source have also been detected by 



Bujarrabal et al.l (119941 ) using the IRAM telescope. In the A 1.3 mm window, we for the first 
time detected the J = 5-4 isotopic transitions of C^^S, C^^S, and ^^CS, among which the 
C^^S (5-4) line emission is relatively faint and its measurement is uncertain. The J = 3-2 
isotopic transitions lie in the A 2 mm region, but are below our detection limit. The two main 
lines have a similar width but are narrower than the isotopic lines, suggesting that the main 
lines are optically thick. We obtained the intensity ratios of /(^^CS 5-4)/ 1 {^^CS 5-4)= 0.09 
and /(C^^S 5-4)//(C^^S 5-4)= 0.13, larger than the corresponding values of 0.03 and 0.07 
in IRC+10216. The trend is the same as that found for the /(^^CO 2-l)/I{^^C0 2-1) ratio. 



3.2.3. C2H 



Fukasaku et al.l (119941 ) detected C2H in this objec t thro ugh the N = 1-0 transition. 
Their detection was further confirmed by IWoods et al.l (120031 ). The N = 3-2 transition is 
split in six hyperfine-structure lines grouped in three fine-structure groups, which lie in our 
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survey range. The two strong main components, 87/2-25/2 and 85/2-23/2, are clearly detected 
in the 2 mm window. The 87/2-25/2 transition is blended with a considerably weaker SiC2 
line. The faint 85/2-25/2 transition is overwhelmed by noise. We obtained the intensity ratio 
of /(C2H 87/2-25/2)//(C2H 85/2-23/2)= 1.7, in good agreement with the values of 1.4 and 1.5 
found in IRC+10216 and CIT6, respectively. 



3.2.4. C4H 



The = 10-9 transition of C4H in this object has been detected by iFukasaku et al. 



(I1994I ). Due to fine structure interaction, every rotational transition of C4H is split into 
two components, which have similar intensities. A total of 12 favorable lines of C4H lie in 
our survey range. The eight C4H lines in the A 2 mm window were clearly detected. The 
four lines in the A 1.8 mm window are too faint to obtain reliable measurements. In Table 
m we only give an uncertain detection of the C4H (24-28 a,b) lines, which is the strongest 
transition of C4H in the A 1.8 mm window. We found that these high- J C4H lines are at 
least 8 times fainter than those in the A 2 mm w indow. In. IRC +10216, however, these C4H 
transitions have a similar strength (see Table 9 of iHe et al.ll2008l ). This cannot be completely 
attributed to different beam dilution effects in the two objects, and can be a reflection of 
the different physical conditions of the two carbon envelopes. 



3.2.5. C-C3H 

The cyclopropynylidyne radical (C-C3H) has not been detected in CRL 8068 before this 
work. Here we report the first tentative detection of this species. Six favorable transitions 
of C-C3H lie in our survey region. Only the strongest two transitions are discovered in our 
spectra with relatively weak strengths. C-C3H has been detected in IRC+10216, whereas 
there is no evidence for the corresponding lines of C-C3H in CIT6. 



3.2.6. CmCN 



There was no report of the detection of CH3CN in this object before this work. IWoods et al. 



(I2OO8I ) estimated the intensity upper limit of the (6i-5i) transition. Here we present the first 
detection of CH3CN in this object. There are a total of 14 favorable CH3CN transitions in 
our survey region. The strongest two fine-structure groups are clearly detected in the A 
1.8 mm window, but the others are below our detection limit. This species was also detected 
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in IRC+10216 and CIT6. 



3.2.1. CN 



The cyanogen radical (CN) is one of the most a bundant molecules in enve lopes of carbon 



rich st ars. CN in CRL 3068 has been discovered by IWoods et al.l (120031 ) and iBachiller et al. 
( 119971 ) in the N = 1-0 and = 2-1 transitions. There are three N = 2-1 fine-structure 



groups of CN present in our survey region. They are split into total of 18 hyperfine structure 
components. All the three groups were apparently observed in the A 2 mm window. The 
weakest one around 226.33 GHz should be a new detection. Our observations suggest the 
intensity ratio of the fine-structure groups around 226.87 GHz an d 226.6 9 GHz to be 1.6, in 
excellent agreement the value of 1.4 obtained by IBachiller et al.l (119971 ). The value is also 
consistent with those found in IRC-l-10216 and CIT6. 



3.2.8. HCN 

The J = 3-2 transitions of HCN and its isotopologues HC^^N and H^'^CN were clearly 
detected with well-defined profiles in the A 2 mm window. The HCN (3-2) transition is the 
second strongest line after the CO (3-2) transition in our survey. A few vibrationally excited 
lines of HCN have been detected in IRC+10216 and GIT 6. The strongest two are detected 
in the spectra of CRL 3068. To our knowledge, this is the first detection of vibrationally 
excited lines in this object. We discovered the isotopologue HC^^N for the first time. The 
intensity ratio of /(HC^^N 3-2)/I(R^^CN 3-2) is 0.027, which is slightly higher than the 
value of 0.016 obtained in IRC+10216. We also obtained the intensity ratio of /(H^'^CN 
3-2)//(H^2QN 3-2)= 0.60, larger than the values of 0.16 and 0.41 found in IRC+10216 and 
GIT 6. The larger isotopologue-to-main line ratios in CRL 3068 are also found for CO and 
CS (see above). Consequently, we infer that the lines in this extreme carbon star is more 
optically thick than those in IRC+10216 and GIT 6. 



3.2.9. HNC 



HNC has been detected in CRL 3068 (ISopka et al.lll989l : lFukasaku et al.lll994J : IWoods et al 



20031 ). This species can be produced through a similar way to HCN. HNC can also be en- 
hanced through ion-molecular reactions or be transferred into HCN in high temperature. 
An abundance comparison of HNC with HCN can provide some insight into circumstellar 
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chemistry. Although no HNC hne hes in our survey region, we conclusively detected its 
isotopologue HN^^C in the J = 3-2 transition for the first time. Our observations suggest 
the intensity ratio /(H^^qn 3-2) / liRN^^C 3-2) to be 19.3, about an order of magnitude 
lower than that in IRC+10216. This implies that HNC has been enhanced in CRL3068. 
HN^^Q ^g^g j^Q^ detected in CIT6. 



3.2.10. HaN 



Cyanoacetylene (HC3N) has been commonly detected in carbor i-rich circumstellar en- 



velopes. The pre s ence of HC.sN in CRL 3068 has been confirmed (I Jewell fc Snvderl Il984 



Sopka et all Il989l : ISuiarrabal et all 1 1994 IPukasaku et~aD Il994j : IWoods et al1l2003h . Nine 
HC3N transitions from J = 15-14 to J = 29-28 lie in the frequency range of our survey. All 
of them were detected in our spectra. Most of these high-J transitions are new detections. 



3.2.11. SiCo 



SiC2 is the species having the most emission lines in our survey. Including possible 
detections and blended features, a total of 23 SiC2 transitions were detected. These SiC2 
lines are relatively weak with peak brightness temperatures less t han 40 mK. W e canri ot find 
previous papers reporting on the detections of SiC2 in CRL 3068. IWoods et al.l (120031 ) failed 
to discover SiC2 in this object and estimated an intensity upper limit of the SiC2 (5o,5-4o,4) 
transition. 



3.2.12. SiO 



SiO in CRL 3068 has been detected by iBujarrabal et al.l (119941 ) in the J = 2-1 and 
J = 3-2 transitions. In our survey the J = 3-2 and J = 6-5 transitions were detected. 
The SiO (6-5) transition was also detected in our spectra of IRC+10216 and CIT6. We 
found that the relative intensity of the SiO transition in CRL 3068 is much lower than 
those in IRC+10216 and CIT6. We obtained the intensity ratio of /(SiO 3-2/ 1 {^'^CO 3-2) 
in CRL 3068 is 0.1, an order of magnitude lower than the values of 1.1 and 1.0 found in 
IRC+10216 and CIT6. Its isotopic transitions ^^SiO and ^°SiO (3-2) lie in the surveyed 
range but are overwhelmed by noise. 



3.2.13. SiS 



The J = 6-5 and J = 5 -4 transitions of SiS in CRL 3068 have previousl y been observed 



by lBujarrabal et all (Il994j ). iFukasaku et all (ll99J), and IWoods et al.l (120031 ). In this survey 
we conclusively detected four SiS transitions with higher J (8-7, 9-8, 13-12, and 14-13) 
transitions. We also report the first possible detection of the isotopologues ^^SiS and ^°SiS. 



3.2.14. Unidentified lines 



Three weak lines detected at 3-8a noise levels remain unidentified. The 264067 MHz 
line probably corresponds to the U264072 line detected in IRC+ 10216 and the other two 
U lines have no corresponding detection in IRC+10216 and CIT6. We also searched the 
NIST frequency table and found that the two features at 245982 MHz and 255108 MHz 
seem to be associated the U lines at 245993 MHz and 255158 MHz detected in Sgr B2(N) 
jNummelin et al.lll998h . 



3.3. Line profiles 



Since the circumstellar envelopes of AGB stars are expanding, molecular line profiles 
can provide significant insight into the nature of the lines as well as the kinematic structure 
of the envelopes. The l ines arising from a spherical envelope might have four characteristic 
profile shapes (see e.g.. lHabing fc Olofssonll2003l : lKwokll2006l ): (a) a fiat-topped line profile 
resulting from optically thin, spatially unresolved emission; (b) a parabolic line profile re- 
sulting from optically thick, spatially unresolved emission; (c) a double-peaked line profile 
resulting from optically thin, spatially resolved emission; (d) a smoothed parabolic (or seem- 
ingly fiat-topped) line profile results from optically thick, spatially resolved emission. Our 
observations show that most of the detected lines are consistent of they being unresolved. 
The profiles of the CO, CS, and HON emission lines appear to be parabolic, indicating that 
they are optically thick. The ^^CO line show a rectangular profile, suggesting that it is likely 
optically thin. Some of the lines, such as vibrationally excited lines of HCN, exhibit fairly 
narrow profiles, probably because they are excited in a compact region close to the central 
star. The different widths for different emission lines seem to suggest that the expansion 
velocity of the envelope is not constant. 
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3.4. Nondetection 

Compared with the spectra of IRC+10216, the nondetected molecular species in our 
survey region include AlCl, AlF, MgNC, NaCN, NaCl, PN, C2S, C3H, C-C3H2, C3N, C3S, l- 
C4H2, H2CO, SiC, and SiN. All the emission lines from these species are weak in IRC+10216 
and are estimated to have an integrated-intensity upper-limit of 0.2Kkms~^ in CRL3068. 
If the line intensities relative to the ^'^CO (2-1) transition in CRL 3068 are the same as those 
in IRC+10216, all these lines are well below our detection limit. Therefore, we cannot draw 
the conclusion that the molecular species in CRL 3068 are less abundant than in IRC+10216. 



4. Column densities and abundances 



For comparison, we have applied the same approach as used in our previous papers 
(IZhang et al.ll2008l . 120091 ) to calculate the column densities and abundances of the molecules 
detected in our spectra. This will minimize the effects of systematic errors when we compare 
the results in different objects. 

The excitation temperatures (Tgx) and column densities (A^) were deduced using the 
standard rotational- diagram analysis. Under the assumption that the lines are optically 
thin, the level populations are in local thermal equilibrium (LTE), and Tex » Ttg, where 
Tbg is the cosmic background radiation temperature (2.7 K), we have the relation. 



m — 

9u 



3k J Tsdv ^ N _ 



(1) 



where Nu, gu, and Eu is the population, degeneracy, and excitation energy of the upper 
level, J Tgdv is the integration of the source brightness temperature over the velocity, S is 
the line strength, /i is the dipole moment, v is the line frequency, Q is the rotational partition 
function. If several transitions arising from levels covering a wide energy range are observed, 
Tex and A^ can be determined using a straight-line fit to Nu/ Qu versus Eu/kTex. Departure 
from the linear relation suggests different excitation mechanisms or misidentification. For 
the calculations, we have corrected the effect of beam dilution through Tg = TjiiOl + 0l)/6l, 
where db (~ 40" and 30" for the ARO 12 m and the SMT respectively) is the antenna full 
beam at half-power, and 9s is the source diameter. Different transitions may have different 9s 
an d it is hard to det ermi ne 9s for each rnolecul e. According to the CO observations reported 
by iNeri et al.l ( 1l998l ) and iTevssier et al.l (120061 ) , we a ssumed 9s = 22" ± 8". The assumed 9s 



Fukasaku et al. 



(I1994I ) and will introduce a ~ 70% uncertainty 



value is similar to that used by 
in the derived column densities. Four molecular species, SiC2, SiS, HC3N, and C4H have 
adequate numbers of detected transitions covering a wide range of excitation energy and the 
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rotational diagrams of these species are given in Figure HI Although in general the transitions 
of these species are well fitted by straight lines, we note that the rotational diagram of SiC2 
indicates a higher Tex for high - J transitions corapared to the lo w- J ones. This is similar to 
what is found in IRC+10216 (lAvery et al.lll992l : iHe et al.ll2008l ). suggesting that excitation 
by infrared radiation may be important for this molecule in AGB envelopes. For the species 
for which the rotational diagrams are unavailable, an average Tex (~ 53 K) was assumed to 
calculate their column densities. 

The molecul ar abundances r espect to molecular hydrogen (/x) were calculated using 
the expression by lOlofssonI (Il996l ). 



fx = 1.7 X 10 



(2) 



where j T^dv is given in Kkms~^, the full half power beam width Oi, is in arc sec, the 
expansion velocity Ve is in kms~^, D is the distance in pc, is the mass loss rate in 
Mq yr~^, Vui the line frequency in GHz, is the statistical weight of the upper level, A^i is the 
Einstein coefficient for the transition, Ei is the energy of the lower level, and Xj e = Ri,e/ (dbD) 
with Ri and Re the inner radius and outer radius of the shell. To aid the reader. Appendix A 
provides full deduction of this formula. For the abundance calcul ation, we first determ ined 
the mass loss rate using the CO (2-1) line. Using the formula of IWinters et al.l (120021 ) and 
assuming fco = 1 x 10~^ and a distance of 1 .05 kpc, we ob t ained ~ 3 x 10~^ yr ~^ 



in good agreement with the recent r esult of ISchoier et al.l (120021 ) and IWoods et al.l (120031 ) 
but a factor of two l ower than tha t of Teyssier et al. ( 20061 ). Whenever available, the Ri and 
Re values given by IWoods et al.l (120031 ) were used. The emission regions of isotopologues 
were assumed to have the same inner and outer radii with those of their main species. The 
Tex values derived from rotation diagram analysis are adopted. For species where more 
than one line is detected, we find that the spreads of the abundances derived from different 
transitions are very small (< 10%), and determine the average values weighted by line 
intensity. Although the adopted Tex value i s significantly larger than the average excitation 
temperature derived by IWoods et al.l (120031 ) for a sample of carbon stars and the excitation 
temperature may not be homogeneous within the envelope, the abundances are essentially 
insensitive to Tex in the case of Ei << kT^^. 

The derived values for Tex, N , and /x are given in Table. [2] . Our calculations are 
based on an optically thin assumption. Therefore, for the optically thick emission, the N 
and /x given in Table [2] should be treated as lower limits. Combining various error sources 
(uncertainties of line intensity, excitation temperature, distance, mass-loss rate, etc.), we 
crudely estimate that the errors of the absolute column densities and abundances amount to a 
factor of ~ 5. Note that the estimated errors are for the case that the basic assumptions of the 
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method, i.e. optically thin emission, are met. If the extent of a given species is significantly 
different from that of CO, its column density may bear a larger error than estimated here. 
It is impossible to accurately determine the extent of each species without interferometric 
observations. In most o f the cases, the sel f-shielding molecule CO is more extended than 
other molecules (see, e.g.. Woods et al. If the source size is overestimated by a factor of 

two, the resultant column density will been underestimated by a factor of four. Accordingly, 
it is possible that the actual errors could rise above a factor of ten. Nevertheless, the column 
density ratios and the abundance ratios should be much more reliable. 



Table [2] also gives a comparison of our results with those derived by iFukasaku et al. 



( I1994J ) and IWoods et al.l (120031 ) . While more molecular species are detected in our survey, 
for the molecules that are commonly detected (but no necessarily from the same transitions). 
Table [2] shows that the and fx values are in agreement within one order of magnitude. 
Such discrepancies are understandable considering different instruments and methods used 
for the observations and abundance analysis. 



Discussion 



5.1. Isotopic species 



Theory predicts that during the AGB phase the nucleosynthesized products in the in- 
ner shell are dredged up to th e stellar sur face, resulting in changing the isotopic ratios 
of circumstellar envelopes 



see 



Busso 2006 



for a recent review). Extensive observations 
have shown that the is otopic composition in carbon-rich stars is remarkably non-solar (e.g. 
Cernicharo et al.ll2000l ). providing substantial tests for nucleosynthesis theory of AGB stars. 
Such non-solar isotopic rati os also make possible the identification of stellar solid-state ma- 
terials in the Solar System (jZinnerlll998l ). 



In Table [3|, we give the isotopic ratios (or their lower limits) of carbon, oxygen, nitrogen, 
sulfur, and silicon in the extreme carbon star, CRL3068. This table also lists the isotopic 
ratios for CIT6, IRC+10216, and the Sun. A comparison of isotopic ratios in these objects 
allow us to study the nucleosynthesis and dredge-up processes in diverse physical conditions. 



5.1.1. Carbon 



The ^^C/^^C isotopic ratio is a good indicator of interior nucleosynthesis processing. The 
dredge-up processes can greatly enhance ^^C in the surface. Moreover, the cool bottom pro- 
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cessing and hot bottom burning can further decrease the C/ C ratio (e.g. lBoothroyd fc Sackmann 



1999 



Frost et aLlll998l ). 



We detected three ^^C-bearing species, ^^CO, ^^CS, and R^^CN. However, the CO/^^CO, 
CS/^^CS, and HCN/H^'^CN ratios only give lower limits of the ^^C/^^C because the main 
lines appear to be optically thick. On the other hand, if the ^^S/^^S isotopic ratio were 
known, we could use the optical thin ^^CS and C^^S lines to obtain the ^^C/^^C ratio. Based 
the fact that the nucleosynthes is in AGE stars i s expe cted not to affect the abundances 
of sulfur isotopologues and that ICernicharo et al.l (I2OOOI ) has found the sulfur isotopic ratio 
in IRC+10216 to be close to solar, we can reasonably assume that the ^^S/'^^S ratio in 
CRL3068 is solar (22.5). From this assumption, we obtained a ^^C/^^C ratio of 29. 7. Using 
a full radiative transfer analysis of ^^CO and ^^CO line emission, IWoods et al.l (120031 ) derived 
the ^^C/^^C ratio of 30, lending strong support to our result. The ^^C/^^C ratio in CRL 3068 
is significantly lower than the solar value, suggesting the presence of extra mixing processes. 

As a consequence of dredge-up processes, the C/0 abundance ratio in circumstellar 
envelopes decreases along with the ^"^C/^^C ratio. As shown in Table [21 the ^^C/^^C ratio in 
CRL 3068 is about a factor of 1.5 lower than those found in CIT 6 and IRC+10216, suggesting 
that CRL 3068 has a higher C/0 ratio. Our observations, therefore, provide an opportunity 
to study circumstellar chemistry in extremely carbon rich environment. 



5.1.2. Nitrogen 



As a star ascends the AGB, the ^^N-rich and ^^N-poor material processed by the CNO 
cycle in the intershell is dredged up to the stel lar surface, resulting a high ^"^N/^^N abundance 



ratio in the envelope. IWannier et al.l (Il99ll ) determined the lower limits of ^^N/^^N in a 



sample of carbon stars and found that in most of the cases they were larger than the terrestrial 
ratio. 

The nitrogen isotopologue, ^^N, has been detected through a faint emission line from 
HC^^N. As its main line is optically thick, the HC^^N/HC^^N abundance ratio gives a lower 
limit of the isotopic ratio of 45. The optically thin species H^'^C^^N and H^^C^^N may also 
be applied to derive the ^''N/^^N ratio. Assuming a ^^C/^^C ratio of 29.7 ( see above), we 
obtained the ^"^N/^^N ratio of 1099, a factor of 4 higher than the solar value. IWannier et al. 



( I1991I ) obtained a lower limit of ^^N/^^N> 500 for CRL 3068, in good agreement with our 
result. The deduced nitrogen isotopic ratio is in accord with the predictions of stellar models. 
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5.1.3. Oxygen 



Nucleosy nthesis mode ls predict enhancement of ^^O in the circumstellar envelopes around 
carbon stars (iBussd l2006l ). Previous observations have shown that the C^^O/C^^O ratios 
i n carbon-rich envelopes are significantly hi gher than the terrestrial and interstellar values 
dWannier fc Sahailll987l : IXahane et aLlll992h . 



The weak and C^^O (2-1) transitions were detected in our survey, allowing us 

to determine the oxygen isotopic ratios in CRL3068. The abundance ratios of CO and its 
isotopologues gives lower limits of isotopic ratios, as can be seen in Table |3l Under the 
assumption of ^^C/^^C= 29.7, we used the optically thin species ^^CO, C^^O, and to 
obtain the ^^O/^^O, ^^O/^^O, and ^^O/^^O ratios of 668, 472, and 0.7, respectively. The 
16q^17q g^j^j 17q^18q radios bear a large uncertainty, but clearly point to an enhancement 
of ^^O. The ^^O/^^O ratio is c omparable with the s olar value. These isotopic ratios are well 
within the ranges reported by iKahane et al.l (119921 ) for a sample of carbon-rich envelopes. 



5.1.4. Sulfur and silicon 



The isotopic compositions of sulfur and silicon are hardly affected by nucleosynthesis of 
AGB stars although they might be slightly changed by neutron capture. This is consistent 
with the observa tions that the S a.nd Si isotopic ratios in IR C-l-10216 and CIT6 are close to 
the solar values (jCernicharo et al.l 120001 : IZhang et al.ll2009l ). 



We detected the S and Si isotopologues ^^CS, ^^CS, ^sSiS, and ^oSiS in CRL3068. The 
main line CS (3-2) is optically thick, while the few SiS transitions are likely to be optical 
thin. Consequently, we derived the ^^S/^'^S , ^^Si/^oSi, ^^Si/^oSi, and ^sSi/^^Si ratios, and the 
lower limit of the ^^S/^^S ratio, as listed in Table O In IRC+102 16 and CIT6, the Si isotopic 



ratios were also obtained from the SiO and SiC2 isotopologues (IHe et al.ll2008l : IZhang et al. 



20091 ). which, however, were not detected in the spectra of CRL3068. We found that these 
S and Si isotopic ratios are compatible with the solar values, and there is no significant 
deviation between the S and Si isotopic ratios in CRL3068 and those in IRC+10216 and 
CIT6. 



5.2. Remarks on the comparison between CRL 3068 and IRC-|-10216 



CRL 3068 and IRC+10216 are the prototy pes of extreme ca rbon stars, which are highly 
evolved heavily obscured by thick dust shells (IVolk et al.lll992l ). Nevertheless, there exist 
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some differences bet ween tfie two object s. CRL 3068 fias a colder IR color temperature than 



IRC+ 10216 (e.g. see lOmont et al.lll993l ). Based on the (ISO) observations. lJusttanont et al. 



(119981 ) have determined the color temperatures of CRL 3068 and IRC+10216 to be of 320 
and 540 K, respectively. The ISO spectra show that the SiC 11.3 /im feature exhibits a 
weak flat-topped emission profile in IRC+10216 while in CRL 3068 it is in absorption. This 
suggests either that CRL 3068 has a thicker dust envelope, or the abundance of the SiC 
grains is higher. Both the color temperature comparison and the SiC feature are consistent 
with the suggestion that CRL 3068 is in an even more evolved state than IRC+10216. 

In Figure [5l we compare the integrated intensity ratios of the lines detected in CRL 3068 
and IRC+10216. The average line ratio between the two objects is 0.06 with a standard 
deviation of 0.07. While there is a large scatter of ratios between species, for individual 
species the ratios are quite consistent. For example, the intensity ratios of SiC2 lines have an 
average of 0.022 with a standard deviation of 0.005. Note that in this a nalysis, we have not 
taken into account the beam dilution effect. Using the same method in IZhang et al.l (120091 ) 
and assuming that CRL 3068 and IRC+10216 have diameters of 22" and 30" respectively, 
we estimate that these intensity ratios would increase by a factor of 1.6 and 1.4 for the ARO 
12 m and the SMT data, respectively if the beam dilution effect is corrected foilfj. CRL 3068 is 
about 8-10 times more distant than IRC+10216, and thus might intrinsically have a stronger 
average line intensity when the distance factor is taken into account. 

Referring to Figure 0, if taking the average line intensity ratio as a reference, one can 
find that the emission from SiO, SiC2, SiS, CS, HCN, and their isotopologues seems to be 
depleted in CRL 3068 compared to those in IRC+10216, while HN^^C, C-C3H, HC3N, and 
CII3CN might be enhanced. The situation for other species is ambiguous. Although the 
classification is not very strict, it reflects that CRL 3068 has different chemical composition, 
excitation conditions, and optical depths with IRC+10216. A detailed discussion of the 
implication for circumstellar chemistry will be given in the following section. 



Our results are in contrast to that presented in lWoods et al.l (120031 ) . who found that for 
a sample of carbin-rich circumstellar envelops, the abundance discrepancies are typically less 
than a factor of five and the two obje c ts, IR C+10216 and CRL 3068, are chemically similar. 
However, in CRL 3068, IWoods et al.l (120031 ) did not detect the species having the largest 
abundance discrepancies found here (e.g. SiO and c-C,sH). Recent stu dies of SiO and SiS i n 
samples of carbon stars have been presented by lSchoier et al.l (120061 ) and lSchoier et al.l (120071 ) . 
who derminated the molecular abundances through a detailed radiative transfer modelling of 



^According to the modelling by I Woods et al.l (|2003f ). CRL 3068 has almost same CO envelope size as 
IRC+10216. If this is the case, the intensity ratios are not affected by the beam dilution effect. 
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the SiO and SiS lines observed by them and some other authors. Their results suggest that 
the fractional abundances of SiO and SiS in carbon stars are in good agreement with each 
other. If the SiO and SiS abundances in CRL3068 are similar to those in IRC+10216, our 
observations would suggest a significant enhancement of carbon chain molecules in CRL 3068, 
as shown in Figure [51 Nevertheless, we cannot completely rule out the possibility that the 
Si-bearing molecules in CRL 3068 are actually depleted. 



5.3. Implication for circumstellar chemistry 

The spectra of CRL 3068 are characterized by the wealth of C-bearing molecules. Among 
the newly detected species, only two are non-C-bearing molecules (SiO and SiS). From Fig- 
ure [5l we can see that CRL 3068 has generally higher intensity ratios of C-bearing molecular 
lines to non-C-bearing ones compared to IRC+10216. This seems to lend support to the 
idea that CRL 3068 is more evolved and more carbon rich than IRC+10216. 

It is intriguing that according to Figure [5] CRL 3068 shows enhancement of HN^^C and 
depletion of H^^CN compared to IRC+10216. Assuming that carbon isotope fractionation 
is insignificant for the two species and the HN^^C/HNC abundance ratio is equal to the 
H^^CN/HCN abundance ratio, we can infer that CRL 3068 has a higher HNC/HCN abun- 
dance ratio than IRC+10216. The HN^^C {3-2)/Il^^CN (3-2) integrated intensit y ratio in 



CRL 3068 is 0.05, about an factor of 13 larger than the value in IRC+10216 (see iHe et al. 



20081 ). HCN and HNC can be produced through a similar way, i.e. dissociative recombination 
of HCNH"*", and can convert to each other at certain physical conditions. At high temper- 
ature, HNC can be reprocessed into HCN through hydrogen exchange reaction. Hence, 
the higher HNC abundance in CRL 3068 might be related to its extreme cold environment. 
Furthermore, HNC can be destructed by atomic oxygen through the reaction 



HNC + O ^ NH + CO. 



(3) 



Therefore, the higher HNC/HCN abundari ce ratio probably reflects the extremely carbon- 
rich nature of CRL 3068. ISchilke et al.l (119921 ) suggest that the HNC/HCN abundance in- 
creases with increasing dust depletion factor because reaction (1) can be suppressed as a 
result of freeze-out of heavy elements onto dust grains. If this is the case, the extreme 
carbon star CRL 3068 may have a larger dust depletion factor than IRC+10216. Our obser- 
vations also suggest that the column den sity ratio of HNC/HCN in CRL 3068 is about one 
order magnitude lo wer than that i n PNs (jJosselin fc Bachilleiil2003l ). This is consistent with 
the observations of iHerpin et al.l (120021 ) who found that HNC is enhanced with respect to 
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HCN at post-AGB stage due to certain chemical processes in the photodissociation region. 



The chemical processes of CN are interconnected with those of HCN and HNC CN 
can be produced through photodissociation of HCN or HNC. On the other hand, CN can 
be reprocessed into HCN or HNC through reaction with molecular hydrogen. Very efficient 



photo dissociation has been found in the highly evolved AGB envelope CIT6 (jZhang et al. 



20091 ). In CRL3068, we do not find evidence showing enhancement of CN. C2H and HC3N 



can be reprocessed by the same chemical precursor, C2H^. The abundance rati o of the two 



species can provide a test for the chemical formation path (jWootten et al.lll980l ). We found 
that the abundance ratio / (C2H )/f(HC.ciN) in CRL 3068 is 7.7, lying within the range of 3-10 
found by IWootten et al. jl980h in a wide variety of interstellar clouds and being consistent 
with gas-phase reaction scheme. 

Among the most interesting finding in this study is the discovery of the cyclic molecule c- 
Cs H in CRL 3068. This follows the detection of another cyclic molecule C3H2 in this source 
by IWannier et al.l (jl99ll ). We further note that HC3N and CH3CN seem to be slightly 
enhanced in CRL 3068. In the proto-planeta ry nebula CRL 618 , cyanopolyynes chains are 
formed through the polymerization of HCN (jPardo et al.ll2005l ). Our observations lead us 
to believe that CRL 3068 may be more evolved than IRC+10216 and cyclic and longer 
chain molecules are quickly being synth esized in its envelope. Since aromatic molecules 
emerge shortly after the end of the AGB (jKwokll2004l ). the possible role that ring and chain 
molecules play in the synthesis of aromatic compounds is an important question in the study 
of circumstellar chemistry. 

The interactio n between stellar winds from AGB stars usually forms shock waves. 
Woods et al.l (120031 ) have found that shocks alter the chemical compositions (e.g. CS and 
SiO) in the inner regio ns of some AGB envelopes. According to a non-equilibrium chemi- 
cal model presented by lWillacy fc Cherchnefil (Il998l ). shocks can also strongly enhance SiO, 
which is a commonly observed species in carbon rich envelopes. However, we do not find 
enhancement of SiO in CRL 3068, probably suggesting that freeze out onto dust grains is 
significant in the extremely carbon rich envelope. 



6. Summary and conclusions 

In this molecular line survey toward the extreme carbon star CRL 3068 in the A 1.3 mm 
and A 2 mm windows, 72 molecular lines from 23 species were detected with three lines re- 
main unidentified. The spectra are dominated by carbon-bearing molecules. The species 
C-C3H, CH3CN, SiCs, C^^O and C^^O, HC^^N, BN^^C, C^^S, C^^S, ^^CS, ^sSiS, and ^oSiS 
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are detected for the first time in this object. From these observations, we have derived the 
chemical abundances and isotopic ratios of the molecular species. A comparison between 
these observations with the spectra of other AGB envelopes obtained previously using the 
same observations settings allows us to rigorously investigate the chemical processes in cir- 
cumstellar envelopes. Compared to the archetypal carbon star IRC+10216, CRL 3068 is 
more carbon rich and shows an enhancement of C-C3H and HN^^C. The overall assessment is 
that CRL 3068 is a more evolved object on the AGB than IRC+10216. However, it should 
bear in mind that our conclusions are mainly based upon the comparison of relative line in- 
tensities which are assumed to directly correspond to the molecular abundances. To obtain 
more reliable absolute abundances and to make stronger statements about the chemical evo- 
lution in circumstellar envelops, interferometric observations and sophisticated treatments 
of radiative transfer processes are required. 

This study is a part of a long-term project of investigating the evolution of circumstel- 
lar chemistry in the late stages of stellar evolution using the A 2 mm and A 1.3 mm spectra 
of a sample of AGB stars, PPNs, and PNs. Such a systematic study of the gas-phase 
molecules presents an important step to understand the synthesis of more complex aromatic 
and aliphatic compounds in evol ved stars, and their role in the distribution of organic com- 



pounds throughout the Galaxy (IKwokl |2009| ) . A detailed study of chemical processes in 



different physical conditions and evolutionary stages will be reported in the future. 
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A. Olofsson formula for molecular abundance 



For the reader's convenience, we deduce the formula presented by lOlofssod (119961 ) for 
the calculation of molecular abundances respect to H2. The basic assumptions include: 1) 
all molecular emission originates from a spherical shell with inner radius and outer radius 
Tg; 2) the expansion velocity Vf.-, the mass loss rate M, and the excitation temperature Tex 
are constant along the radius 3) the molecular abundance respect to H2 is homogeneous in 
the envelope. 

For a Gaussian beam measuring a molecular line, the main beam temperature at a given 
velocity is obtained by 

b 

where p is the angular distance to the center, 0^ is the half-power beam width, and Ts{v,p) 
is the source brightness temperature at v and p. For a given volume element at a distance 
to the center of the envelope r, we have 



p=L^l-{v/vef, (A2) 

where D is the distance between the antenna and the source and v is the velocity along the 
line of sight. From the radiative transfer equation, the source brightnesss temperature is 
given by 

T.(t;,p) = rex[l-e--(^'P)], (A3) 
where t{v,p) is the optical depth. Substituting Equ. A2 and Equ. A3 into Equ. Al, we have 

Tr{v) = 8 In 2Te. r e-4i°2[i-(.K)2].-^ _ [v/vef] [1 - e-"("'")]a;da;, (A4) 

J Xi 

where we have defined x = r/{6hD) and thus Xi^e = fi^e/{ObD). 
Accroding to radiative transfer equation, t{v^x) is given by 

_ ^-hu/kT.^y-EUkT..^^^ (A5) 



c^guAuin{x) 



where Qu is the weight of the upper level, A^i the Einstein coefficient, h the Planck constant, 
k the Boltzman constant, u the line frequency, n{x) the number density of the species, QiT^^) 
the partition function, Ei the energy of lower level, and Sz the path length element along 
the line of sight. From the relation z = r{v/ve) = Dp[l — {v / VeY]~^'^ {v / Ve) , we have 

Sz = ^[1 - {v/ve)T^-''6v = ^[1 - {v/vey]-^6v. (A6) 
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Prom Equs, A5 and A6, we obtain 



where 



^i^^^)-^[^-i^/^efr'r{0,l), (A7) 
The number density for a molecular species is given by 

= 2^ 

or 

where /x is the abundance relative to H2, m-n^ is the mass of a hydrogen molecule, and 

is the mass loss rate for H2. Substituting Equs. A7-A9 into A4 and assuming t{v,x) « 1, 

we obtain the molecular abundance relative to H2 
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Fig. 1. — Spectrum of CRL3068 in the frequency range 130-163 GHz obtained with the 
ARO 12 m telescope. The spectra have been smoothed to a frequency resolution of IMHz. 
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Fig. 1. — continued. 
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Fig. 1. — continued. 
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Fig. 1. — continued. 
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Fig. 1. — continued. 
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Fig. 2. — Spectrum of CRL3068 in the frequency ranges 219.5-267.5 GHz obtained with the 
SMT 10 m telescope. The spectra have been smoothed to a frequency resolution of 3 MHz. 
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Fig. 2. — continued. 




Fig. 2. — continued. 
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Fig. 3. — Cumulative number o f observed lines exce eding a give n valu e. The spectra of 
IRC+10216 and CIT 6 are from Ifle et aP J2008h and Izhang et all J2009h . respectively. All 
the spectra were obtained using the same telescope settings in the same frequency ranges. If 
the distance to an object is increased by a factor of two, the curve would move toward left 
by a length denoted by the shaft of the arrow in the lower left-hand corner. 
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Fig. 4. — Rotational diagrams for the detected species in CRL 3068. 
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Fig. 5. — Integrated intensity ratios of the lines detected in CRL3068 and those detected 
in IRC+10216. The dotted hne represents the average value. 
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Table 1. Molecular transitions detected in CRL3068. 



Species 


Transition 


Frequency" 


rms 






AVpwHM 


(MHz) 




(upper-lower) 


(mK) 


(K) 


(K km/s) 


(km/s) 


CO 


J=2-l 


230538.0 


4 


1.893 


35.942 


17.7 




J=2-l 


220398.7 


4 


0.221 


5.149 


21.2 


C O 


J=2— 1 


224714.4 


4 


0.009 


0.23: 






J=2— 1 


219560.4 


4 


0.016 


0.313 


25.0 


cs 


J=3-2 


146969.0 


11 


0.143 


2.658 


18.4 




J=5-4 


244935.6 


4 


0.201 


3.821 


17.3 


O b 


J =5— 4 


A on 1 Q 

z4zyio.o 


3 


0.010 


0.17; 






J =5— 4 


Z41Ulb. 1 


3 


0.025 


0.485 


21.2 


Ob 


J=5-4 


231221.0 


4 


0.013 


0.325 


33.5 




IN=d— 2 J=7/2— 5/2 


262005. 3* 


7 


0.130 


3.233 


22.3 




JN=3— 2 j=o/ 2—6/ 2 


262066.1 


7 


0.079 


1.950 


31.6 




N=14— 13 a 


100010 T 

133213.7 


14 


0.052 


1.274 






JN=14— lo b 


1 ooo cro 1 

loozoz.l 


14 


0.035 


1.100 






iN=lo— 14 a 


14Z / Zo.o 


11 


0.030 


0.452 


27.2 




"NT 1 c; 1 /I K 


142 m (^.o 


11 


0.029 


0.553 


27.3 




IN — io— lo a 


10ZZ4o.O 


y 


n nop. 


U.4oy 


2o.o 






1 c;onon i 


9 


0.029 


0.508 


27.7 




N=17— 16 a 


1 1 f7TZ 1 

lOlYOO.l 


8 


0.036 


0.715 


20.8 




iN = i/— lb b 


Ibl (>)0.b 


8 


0.041 


0.716 


22.2 




"XT /I Q 1^ 


000 Q /I C 

2zoo4o.O 


4 


0.010 


0.26: 






M /I Q K 

iN=z4— 2v5 b 


000 Q n 
225oo( .U 


4 


0.009 


0.32: 






5(l,4j-4(l,3j 11/2-9/2 


2o2by7.o 


4 


0.010 


0.23: 






5(1,4)-4(1,3) 9/2-8/2 


252801.0 


4 


0.010 


0.31: 






Jk=12i-11i 


22U74o.U 


4 


0.014 


0.291 


13.2 




Jk=12o-11o 


22U747.2^ 


4 










Jk=142-132 


257507.6 


4 


0.015 


0.264 


26.0 




Jk=14i-13i 


257522.4 


4 










Jk=14o-13o 


257527.4^ 


4 










iN(^J,r — 2^^0/2, i/2_)—H^o/2,i/2_) 


zzozo 1 .-y: 




U.UiU 


U.2U. 






N(J,F)=2(3/2,1/2)-1(3/2,3/2) 


226298.9* 


4 










N(J,F)=2(3/2,3/2)-l(3/2,l/2) 


226303.0* 


4 










N(J,F)=2(3/2,3/2)-1(3/2,3/2) 


226314.5* 


4 










N(J,F)=2(3/2,3/2)-l(3/2,5/2) 


226332.5* 


4 


0.015 


0.594 


22.5 




N(J,F)=2(3/2,5/2)-1(3/2,3/2) 


226341.9* 


4 










N(J,F)=2(3/2,5/2)-l{3/2,5/2) 


226359.9* 


4 










N(J,F)=2(3/2,1/2)-1(1/2,3/2) 


226616.6* 


4 


0.016 


0.405 


23.8 




N(J,F)=2(3/2,3/2)-l(l/2,3/2) 


226632.2* 


4 










N(J,F)=2(3/2,5/2)-l(l/2,3/2) 


226659.6* 


4 


0.036 


1.049 


21.8 




N(J,F)=2(3/2,l/2)-l(l/2,l/2) 


226663.7* 


4 










N(J,F)=2(3/2,3/2)-l(l/2,l/2) 


226679.3* 


4 










N(J,F)=2(5/2,5/2)-l(3/2,3/2) 


226874.2* 


4 


0.088 


2.399 


21.9 




N(J,F)=2(5/2,7/2)-l{3/2,5/2) 


226874.8* 


4 










N(J,F)=2(5/2,3/2)-l(3/2,l/2) 


226875.9* 


4 










N(J,F)=2(5/2,3/2)-l(3/2,3/2) 


226887.4* 


4 
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Table 1 — Continued 



Species 


Transition 


Frequency" 


rms 


Tr 




AVfwhm 


(MHz) 




(upper-lower) 


(mK) 


(K) 


(K km/s) 


(km/s) 




N(,J,F)=2(5/2,5/2)-l (3/2,5/2) 


226892.1* 


4 










N(J,F)=2(5/2,3/2)-l(3/2,5/2) 


226905.4* 


4 








HCN 


J =3-2 


265886.4 


7 


0.831 


14.315 


15.8 




V2=l^'^ J=3-2 


265852.8 


7 


0.034 


0.401 


11.5 




1^2=2° J=3-2 


267243.2 


7 


0.036 


0.231 


14.5 


HO IN 


J=3-2 


258154.7 


4 


0.014 


0.233 


22.3 


TT 1 3 /^M 


J =3-2 


259011.8 


5 


0.495 


8.582 


15.8 


T T "V T 1 ^ 


J =3-2 


261263.4 


7 


0.021 


0.445 


23.7 


HC3N 


J=15-14 


136464.4 


11 


0.082 


1.787 


21.4 




J=16-15 


145560.9 


10 


0.053 


0.953 


21.4 




J=17— 16 


154657.3 


9 


0.020 


0.45: 






J=18-17 


163753.4 


10 


0.075 


1.611 


21.8 




J=25— 24 


z27418.9 


4 


0.027 


0.530 


20.5 




J=26-25 


236512.8 


6 


0.031 


0.672 


25.3 




J=27-26 


245606.3 


4 


0.031 


0.541 


21.3 




J=28— 27 


254699.5 


5 


0.021 


0.485 


19.3 




J=29— 28 


263792.3 


6 


0.021 


0.310 


26; 


SiC2 


6o,6-5o,B 


137180.8 


10 


0.022 


0.424 


21: 




62,5-52,4 


14U92U.1 


9 


0.025 


0.483 


20: 




62,4-52,3 


145325.8 


9 


0.029 


0.560 


20.6 




7o,7-6o,6 


158499.2 


7 


0.033 


0.617 


25.0 




10o,io-9o,9 


220773.7 


4 


0.019 


0.380 


18.3 




92,7-82,6 


222009.4 


5 


0.022 


0.400 


29.0 




102,9-82,8 


232534.1 


4 


0.018 


0.411 


29.7 




108,2-98,1 


234534.0^ 


6 


0.012 


0.222 


29.3 




108,3-98,2 


234534.0 


6 










106,5-96,4 


235713.0 


5 


0.024 


0.415 


15: 




106,4-96,3 


235713.1 


5 










104,7-94,6 


237150.0 


5 


0.019 


0.417 


30: 




104,6-94,5 


237331.3 


5 


0.024 


0.35: 


15: 




llo,ii-10o,io 


241367.7 


3 


0.020 


0.439 


18.6 




102,8-92,7 


247529.1 


5 


0.022 


0.544 


22.8 




112,10-102,9 


254981.5 


5 


0.030 


0.466 


28: 




118,3-108,2 


258065.0* 


4 


0.014 


0.189 


16.2 




118,4-108,3 


258065.0* 


4 










116,6-106,5 


259433.3* 


5 


0.019 


0.502 


23.0 




116,5-106,4 


259433.3* 


5 










114,8-104,7 


261150.7 


7 


0.016 


0.33: 






114,7-104,6 


261509.3 


7 


0.015 


0.28: 


22.4 




12o,12-llo,ll 


261990.7$ 


7 


0.130 


3.233 


22.3 


SiO 


J=3-2 


130268.6 


12 


0.026 


0.976 


22.1 




J=6-5 


260518.0 


7 


0.034 


0.537 


15.7 


SiS 


J=8-7 


145227.0 


9 


0.078 


1.517 


20.3 




J=9-8 


163376.7 


10 


0.119 


2.127 


19.5 
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Table 1 — Continued 



Species 


Transition 


Frequency" 


rms 


Tr 


jTndv'' 


AVfwhm 


(MHz) 




(upper-lower) 


(mK) 


(K) 


(K km/s) 


(km/s) 




J=13-12 


235961.1 


6 


0.101 


2.188 


18.8 




J=14-13 


254102.9 


5 


0.121 


2.265 


18.4 


29SiS 


J=13-12 


231626.7 


4 


0.013 


0.243 


30: 


30SiS 


J=9-8 


157579.7 


9 


0.025 


0.44: 






J=15-14 


262585.0 


7 


0.021 


0.101 


19.4 


U 




245982.0 


5 


0.017 


0.605 


14.7 






255108.0 


5 


0.041 


0.817 


16.0 






264067.0 


5 


0.016 


0.175 


13.9 



"♦-unsolved hyperfine structure lines; $ blended with other species. 

'The colon indicates uncertain detections. The pointing uncertainties induced errors 
are up to 20%-30%. 
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Table 2. Excitation temperatures, column densities and abundances respect to H2". 



Species 


Tex (K)*- 


N (cm 




/x 




This paper'^ 


¥94" 


This paper'^ 


■^¥03" 


CO 




6.20(16) 




1.5(-4) 








1 17^'1^^^ 

1. 1- 1 \^10J 




Z.l(-C)) 








5.58(14): 




1.2(-6): 








7.66(14) 




1.7(-6) 




cs 




3.24(13) 




5.8(-8) 


3.7(-7) 






1.38(12): 




2.7(-9): 








3.95(12) 




7.8(-9) 








2.91(12) 




5.9(-9) 




C2H 




1.82(14) 


1.5 ± 1.1(14) 


2.7(-7) 


5.7(-6) 


C4H 


38 ±18 


3.52(14) 


7.0(14) 


6.5(-7) 




C-C3H 




1.78(13): 




3.6(-8): 




CH3CN 




5.77(12) 




8.1(-8) 


< 1.2(-7) 


CN 




1.93(14) 




3.9(-7) 


4.6(-7) 


HCN 




4.17(13) 


6.0(14) 


4.5(-8) 


6.3(-6) 






5.49(11) 




1.0(-9) 








2.01(13) 




3.7(-8) 


3.0(-8) 


HNi^C 




1.26(12) 




2.3(-9) 


< 1.4(-8) 


HC3N 


46 ± 15 


1.78(13) 


1.9 ±1.7(14) 


3.5(-8) 


5.0(-7) 


SiC2 


76 ±36 


3.49(13) 




3.8(-8) 


< 7.5(8) 


SiO 




1.64(12) 




3.0(-9) 




SiS 


52 ± 16 


5.69(13) 


5.4 ±9.0(13) 


1.5(-7) 


6.7(-7) 


29SiS 




6.38(12) 




1.3(-8) 




30SiS 




2.83(12) 




5.2(-9) 





'^x{y) represents x x lO^*; 

* A constant excitation temperature of 53 K was assumed for the species for which the 

rotation diagrams cannot be obtained; 

"^For the species with optically thick emission (e.g. CO and HCN), this gives the 
lower limits; 

'*F94: from Fukasaku et al. (1994); W03: form Woods et al. (2003). 
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Table 3. Isotopic abundance ratios. 



Isotopic ratio 


CRL 3068 
Sepcies 


Value 


CIT6« 


IRC+10216'' 


Solar^ 




19/^34n /13/^39n 

C S/ s 


29.7'^ 


45.4 


45 


89 






5.6"^ 










1 9/^o /IS /^o 

i^og/idog 


9.8^ 










H^^CN/H^CN 


1.2"= 










Jjl3Ql4j^/J|12Ql5p^ 


1099-^ 






272 






45'' 








160/170 


130160/120170 


668:^ 


890 


967 


2680 






125:^ 








160/180 


130160/120180 


472-'^ 




1172 


499 




C^^o/Ci^O 


88^ 








170/180 




0.7: 




1.14 


0.2 


32g/34g 




7.4^ 


6.7^ 


18.9 


22.5 


33g/34g 


C33S/C34S 


0.3: 


0.2: 


0.19 


0.18 


29Si/30Si 


29SiS/30SiS 


2.5 


1.0 


1.46 


1.52 


28gi/30gi 


28SiS/30SiS 


28.8 


8.8^ 


24.7 


29.9 


28Si/29Si 


28SiS/29SiS 


11.5 


8.9^ 


17.2 


19.6 



''From Zhang et al. (2009); 

*From He et al. (2008) except the ^^C/^^C ratio which was taken from Cer- 
nicharo et al. (2000); 

''From Lodders (2003); 

'^Assume that the ^''S/^^S ratio is solar; 

''Should be treated as lower limits due to opacity effect; 

^Adopted: 120/130=29.7. 



